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Abstract 

At low temperature, spheres with a very short-ranged attraction exist as a near-close- 
packed solid coexisting with an almost infinitely dilute gas. We find that the ratio of the 
interfacial tension between these two phases to the thermal energy diverges as the range of 
the attraction tends to zero. The large tensions when the interparticle attractions are short- 
ranged may be why globular proteins only crystallise over a narrow range of conditions. 

PACS: 68.10.Cr, 68.35.-p, 82.70.Dd 



1 Introduction 

The interfacial tension is a useful quantity to know, it not only defines the cost of an interface 
but is also a central feature in the classical nucleation theory of first order phase transitions Q . 
However, the interfacial tension between solid and fluid phases is unknown for all but a very few 
off-lattice microscopic models. Even for the very simple hard-sphere potential a consensus on 
its value has only recently been reached This is despite the fact that we have known the 

bulk phase diagram of hard spheres for thirty years Q . The reason for the lack of calculations of 
surface tensions is due to their difficulty: calculation of the interfacial tension of hard spheres is 
a formidable problem in density functional theory 0, |3| . Here we calculate the interfacial tension 
of hard spheres with a very short-ranged attraction. The limiting case when the range tends 
to zero is the (3q model of Stell || ^| . The phase behaviour of the (3$ model is straightforward, 
if a little peculiar. Above a certain temperature T co u the behaviour is identical to that of 
hard spheres (below close packing |7]]) and below this temperature they phase separate into an 
infinitely dilute gas coexisting with a close-packed solid |5]-||] . 

The bulk solid phase of spheres with a very short-ranged attraction can be described, due to 



its very high density, accurately and simply using a cell theory [1C]. We will extend our previous 
cell theory treatment of the bulk || to the interface in order to calculate the interfacial tension 
analytically. Of course, the interfacial tension between a solid and another phase depends on 
the orientation of the interface with respect to the lattice of the solid. We calculate it when the 
surface of the solid is one of a few of the low index lattice planes. 

The Po limit with its zero ranged attraction is a purely mathematical limit, unobtainable in 
an experiment. However, there are two types of colloidal systems where the attraction has a 
range which is small in comparison to that of the hard repulsive interaction. These are mixtures 
of colloidal particles with either smaller particles |pr|l3 (which may be surfactant micelles [O] ) 



or small polymer coils [16|, and globular protein molecules []T7|-pO||. 



2 Model and bulk phase behaviour 



First, we define the well-known square-well potential. It is the spherically symmetric pair po- 
tential u(r) defined by 

{oo r < a 

-e a<r<a(l + 5) , (1) 

r>a(l + 5) 

where a is the hard-sphere diameter, and r is the separation between the centres of the spheres. 
Here we will always be considering short-ranged attractions, 5 <C 1. The first person to consider 
very short-ranged attractions was Baxter [21] who considered a potential with zero range, 5 = 0, 



and with a well-depth e/T adjusted so that the second virial coefficient was of order unity. 
This model is often termed the sticky-sphere model. Within it the second virial coefficient is 
used as a temperature like variable. However, Stell || showed that the sticky-sphere model 
was pathological, its fluid phase was unstable at all non-zero densities. Therefore, we will not 
consider this model but instead will follow Stell when we take the limit 5 — > 0, thus obtaining 
his Po model ||. 

The bulk phase behaviour of the /3q model is described in Refs. 0-0,0] • If the two limits 
5 — ► and T/e — > are taken such that T > T co n then the equilibrium phase behaviour is 
identical to that of hard spheres. If T < T co u then the behaviour is radically different: a close- 
packed solid coexists with a fluid phase of zero density. This close-packed solid may be either 
face-centred cubic or hexagonal close-packed; both have the same number of nearest neighbours 
and the same maximum density and so will have very similar free energies. Because of this we 
do not specify which one is actually formed. The temperature T co u is [§|,||,|9| 

e ~ ln(l/$) ' 1 ' 

Note that we have used energy units for the temperature, i.e., units in which Boltzmann's 
constant equals unity. 



3 Interfacial tension 

Above T co n the interface between the coexisting fluid and solid phases is identical to that of 



hard spheres; see Refs. |2|,|3, 22 for work on the interface between solid and fluid phases of hard 
spheres. Below T co u the gas-solid interface will be very different. It will be very narrow. The 
free energies per particle (/T) in the coexisting solid and gas phases approach — oo as the (5q 
limit is taken, while at intermediate densities, all densities which are non-zero and below close 
packing, the free energy is much higher 0]. Thus we expect the density will remain almost at 
its value in the bulk solid even in the outermost layer of the solid before dropping abruptly to 
zero. We will calculate the surface tension on the basis that only the outermost layer of the 
solid differs from the bulk. The solid layer beneath it, and the gas phase right up to the solid, 
are assumed to be identical to the bulk solid and gas phases, respectively. 



The free energy of a bulk solid phase can be estimated using a cell theory [10|. A cell theory 
starts from the 1-particle partition function q\ of a particle trapped in a cell formed from its 
neighbouring particles fixed at the positions they occupy in a ideal lattice. The free energy per 
particle a s is then obtained from 

-^=-\nqx, (3) 

where q\ is defined to be in units of a 3 , and a term which is the logarithm of the thermal volume 
of a sphere divided by <r 3 is neglected. For a close-packed solid q\ is, see Ref. M, 



</i.^v>xpk; ( /7) p > (1+ P J /2)3 - (4) 
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p cp = V2a~ 3 is the close-packed density of a close-packed solid of hard spheres. The restriction 
on the minimum density of the solid ensures that the sphere is close enough to all twelve of its 
neighbours to interact via the attraction. Inserting Eq. @ into Eq. (^) yields 

t=- 31u6 -t P> (TTJJW () 

The first term in Eq. (||) is the logarithm of the volume available to the centre of a sphere 
and the second term is (half) its energy of interaction with its twelve neighbours. The volume 
available to the centre of a sphere is not precisely 5 3 but is c5 3 , where c is a prefactor of order 
unity. We have neglected the lnc term in Eq. (||) as it is of order unity whilst the other terms 
diverge in the (3q limit. 

Equation (|5|) gives the free energy per sphere in the interior of the solid phase. The free 
energy per sphere in the outermost layer will be different. It will be given by an expression of 
the form of Eq. (||) but in which q\ is replaced by the partition function of a particle in the 
outermost layer, q\. A particle in the outermost layer has fewer neighbours than in the bulk; 
recall that the coexisting gas is at very low density so there are almost no interactions between 
the outermost layer and the gas. How many fewer neighbours depends on which lattice plane 
forms the outermost layer. We denote the number of missing neighbours by z m , it is equal to 
three for an interface in the 111 plane of a face-centred-cubic lattice. So, the energy of a particle 
in the outermost layer is (6 — z m /2)e. All we now need is the available volume for a particle in 
the outermost layer. The particle can only explore a volume greater than ~ <5 3 at a cost of no 
longer being within 5 of all its remaining 12 — z m neighbours and so increasing the energy. It 
is easy to see that if the particle moves over a distance much larger than 5 in any one direction 
then it can only remain within 5 of two particles: its motion consists of rolling over the surfaces 
of a pair of adjacent spheres and so is restricted to a volume of order 5 2 a. Thus, the entropy 
gain (xT) is — Tln<5 but the energy cost is (10 — z m )e. For the solid to be stable the temperature 
must be below T co u, Eq. (||), and so —Tin 6 < 2e. The entropy gain is then is only greater than 
the energy cost when z m > 8. Similar arguments apply for allowing a particle to move over a 
distance much larger than 5 in two or three directions. This is only favourable when z m > 7 
and z m > 6, respectively. For any flat outer layer z m will be less than six, and the outermost 
particles will rattle inside a volume of order 5 s as they do in the bulk. This is consistent with 
our assumption that only the outermost layer of the solid differs from the bulk. The particles in 
the layer below the outermost layer interact with twelve neighbours and so have the same free 
energy as in the bulk, apart from corrections of order T. 

So, for qf we have 

g|~5 3 exp[(6-W2)e/T] p> {1 ^ /2)3 , (6) 

which gives a free energy difference per particle between the outermost layer and the interior of 
the solid of 

T ~ ~2T 9 (1 + 5/2)3- [ ) 

This difference can be converted into a surface tension by dividing by the area per sphere in the 
outermost layer. (In doing so we are implicitly fixing the surface of tension to be that which 
fixes the surface excess number of particles to be zero [23|.) For example, for a 111 surface the 
area per sphere is ^/3/2a 3 and z m = 3, so 

7m = ^3ecj- 2 ~ 1.73«r~ 2 T < T coU . (8) 

Similarly, for the 110 and 100 surfaces of a face-centred-cubic lattice the areas per sphere are 
^/2cr 2 and a 2 , respectively and the z m 's are 5 and 4. Thus, 7110 = (5/(2^/2))eo"~ 2 ~ 1.77eer~ 2 
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and 7100 = 2ea~ 2 . The 111 surface has the lowest surface tension because it has the lowest ratio 
of number of missing bonds per surface sphere to area per sphere. The surface tensions are all 
of order ea~ 2 . In the (3q limit the ratio e/T is infinite, recall that we are below T co u. Thus, 
the ratio of the surface tension (expressed using the sphere diameter as a unit of length) to the 
thermal energy is infinite. As the range of attraction becomes very small the interfacial tension 
becomes very large. 

The assumptions which underlie the derivation of Eq. (||) should be valid whenever e/T S> 1, 
8 < 0.1 and the fluid phase is highly dilute. The first restriction ensures that the interface is only 
one layer thick, the second that Eqs. (||) and (|6|) are valid, and the third that the outer layer of 
solid does not interact to a significant degree with the fluid phase. Under these conditions, Eq. 
(|8j) is a good approximation for the interfacial tension and this tension is high. 

An interfacial tension of order ea~ 2 is not a surprise. It is what we would obtain if we just 
approximated the surface tension by the energy per unit area needed to pull a block of solid 



apart to create two new surfaces; see the book of Israelachvili 24] where he estimates interfacial 
tensions using just such an approximation. Note that Eq. @ has no explicit dependence on the 
range of the attraction 5, and so is a rough estimate of the low temperature interfacial tension 
even of the Lennard-Jones potential. Then e would be the well-depth of the Lennard-Jones 
potential. 

We now compare our results with the earlier work of Marr and Gast |25,^6|. This work 



was conducted within the Percus-Yevick (PY) approximation for sticky spheres [21|. PY for the 
sticky-sphere model was shown by Stell f|,|6],|| to yield qualitatively incorrect results; it predicts 
vapour-liquid equilibrium at a temperature of e/(ln(l/5)). This is below T co u, Eq. (||), and so 
in fact the fluid phase is unstable at all non-zero densities. Thus the results of Marr and Gast 
| 25| , 26] are for the interface between phases which do not exist. 



4 Discussion and Consequences 

We have calculated the low temperature interfacial tension of hard spheres with a very short- 
ranged attraction and found that its ratio to the thermal energy per unit area is very large. Low 
temperature means below T co u, Eq. (0). In the limit that the range of the attraction tends to 
zero, 5 — > 0, the 0o limit J|], then this ratio diverges. In this limit above T co u the attractive part 
of the interaction has a negligible effect: all the equilibrium properties, including the interfacial 
tension are identical to those of hard spheres. At T co u, the coexisting fluid and solid densities 
change discontinuously to zero and close-packing (p cp ), respectively and the interfacial 

tension jumps to that given by Eq. (^). If the range 5 of the attraction is very small but 
finite then the discontinuity at T co u becomes a narrow temperature range [^?J over which the 
fluid and solid densities at coexistence rapidly decrease and increase, respectively. Below this 
temperature range the density of the fluid is very low and our Eq. (|8|) for the interfacial tension 
will be accurate. 

When the range of the attraction is very small there is only a narrow temperature range 
[27] separating a high-temperature regime in which the spheres are almost hard spheres and 



a low-temperature regime where they are almost at the low-temperature limit. By almost at 
the low-temperature limit we mean that there is a very dilute gas coexisting with a solid with 
a density near close packing, and the interfacial tension between the two phases is then very 
high. This may explain the finding of George and Wilson [28] that there is only a narrow slot in 
effective temperature within which globular proteins can be made to crystallise [p^-20|. (Other, 
not necessarily contradictory, explanations have been proposed by Poon [29], and by ten Wolde 
and Frenkel p6| .) The narrow slot may correspond to the narrow temperature range where 
the coexisting densities and interfacial tension are changing rapidly. Above this temperature 
range the spheres are hard-sphere-like and so only crystallise at high density, above a volume 
fraction of 0.49 — the density at which a fluid of hard spheres coexists with the solid [fl. At 
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this density the dynamics of crystallisation may be slow |3l],|33l due *° a nearby glass transition. 
Below this temperature range the interfacial tension is very large. The free energy barrier to 
nucleation varies as the cube of the interfacial tension, within classical nucleation theory Q. 
Thus, the rate of homogeneous nucleation varies as exp(— 7 3 ) and so is extremely small when 
the interfacial tension is large. We conclude that spheres with a short-ranged attraction only 
crystallise easily from a dilute solution over a narrow temperature range around T co u: above it 
the spheres crystallise only at high density; below it the interfacial tension is large and hence 
homogeneous nucleation is extremely slow. 
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